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O ■ We present a detailed numerical study of the dynamics and evolution of ul- 

trarelativistic magnetohydrodynamic jets in the black hole-disk system under ex- 
treme magnetization conditions. We find that Lorentz factors of up to 3000 are 
achieved and derived a modified Michel scaling (F ~ a) which allows for a wide 
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■ variation in the flow Lorentz factor. Pending contamination induced by mass- 

entrainment, the linear Michel scaling links modulations in the ultrarelativistic 
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wind to variations in mass accretion in the disk for a given magnetization. The jet 
is asymptotically dominated by the toroidal magnetic field allowing for efficient 
collimation. We discuss our solutions (jets) in the context of Gamma ray bursts 
and describe the relevant features such as the high variability in the Lorentz 
factor and how high collimation angles (~ 0° — 5°), or cylindrical jets, can be 
achieved. We isolate a jet instability mechanism we refer to as the "bottle-neck" 
instability which essentially relies on a high magnetization and a recoUimation 
of the magnetic flux surfaces. The instability occurs at large radii where any 
dissipation of the magnetic energy into radiation would in principle result in an 
optically thin emission. 

Subject headings: gamma rays:bursts — magnetic flelds — MHD — ISM:jets and 
outflows 
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1. Introduction 

It is widely accepted that the most conventional interpretation of the observed GRBs 
result from the conversion of the kinetic energy of ultra-relativistic particles (wind) to ra- 
diation in an optically thin region. The particles being accelerated by a fireball mechanism 
taking place near the central engine (Goodman 1986; Shemi & Piran 1990; Paczyiiski 1990). 
The prompt 7-ray emission is probably induced by internal shocks within the wind while the 
afterglow results from the wind external shock interacting with the surrounding medium. 
The Lorentz factor F of the relativistic wind must reach high values (F 10^ — 10^) both 
to produce 7-rays and to avoid photon-photon annihilation along the line of sight, whose 
signature is not observed in the spectra of GRBs (Goodman 1986). The crucial point is 
that for a range of plausible parameters, the prompt occurs above the Compton photosphere 
(Meszaros 2002). 

Whether produced locally or originating from the source, large magnetic fields are re- 
quired to account for the synchrotron and/or inverse Compton emission from a non-thermal 
population of accelerated electrons produced behind the (internal and external) shock waves. 
In the case of the afterglow emission, the magnetic field has to be locally generated by mi- 
croscopic processes (Meszaros & Rees 1993; Wijers, Rees, & Meszaros 1997; Thompson & 
Madau 2000). In the case of the prompt 7-ray emission, such a locally generated magnetic 
field is also usually invoked Rees & Meszaros (1994); Papathanassiou & Meszaros (1996); 
Sari & Piran (1997). Although a large scale field originating from the hidden source could 
play the same role (Meszaros & Rees 1993, 1997; Tavani 1996). 

Support for a high coUimation of these winds is derived from the break in the light curves 
of afterglows for long duration gamma ray bursts (e.g. Stanek et al. (1999)). The brightness 
of the optical transient associated to GRB990123 showed a break (Kulkarni et al. 1999), and 
a steepening from a power law in time proportional to ultimately approaching a slope 

^-2.5 (Castro-Tirado et al. 1999). The achromatic steepening of the optical light curve and 
early radio flux decay of GRB 990510 are inconsistent with simple spherical expansion, and 
well flt by jet evolution. 

The high Lorentz factors invoked, the reference to a magnetic field and the high coUima- 
tion are ingredients suggestive of magnetized jets at play in these phenomena. The extreme 
magnetic fields led to models using a highly magnetized millisecond pulsar. However, due 
to the baryon loading, GRB jets are different from the essentially baryon-free pulsar wind 
case (e.g. Usov (1994)). Magnetohydrodynamic (MHD) jets emanating from the surface of 
an accretion disk surrounding a central object are naturally loaded with baryons and offer 
a more interesting case for GRBs. Necessary currents are maintained given the standard 
baryon/mass loading which guarantees the validity of the MHD approximation out to large 
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distances (up to 10^° cm; Spruit et al. (2001) and Appendix A in this paper). The MHD 
approximation simpUfies the mathematical treatment of the problem but on the other hand 
internal magnetic energy dissipation mechanisms are to be sought that would lead to the 
proper (GRB-like) emission. The MHD approximation breaks down when the flow reaches 
the critical radius beyond which current can no longer be sustained; wc will refer to this 
radius as the MHD radius. Note that for the pulsar wind case, the MHD approximation 
breaks down much earlier, and plasma theories of large amplitude electromagnetic waves are 
applied to explain the prompt emission (Lyutikov & Blackman 2001). Here, we focus on 
MHD jets emanating from the disk-black hole system (Thompson 1994; Meszaros & Rees 
1997). We are interested in solutions that account for the crucial features inherent to GRB 
models such as i) the extreme and variable Lorcntz factors as in the internal shock model, 
ii) the high collimation {9 1° — 10°), iii) the dissipation mechanisms. We start in Sec. 2 by 
introducing the reader to the basic concepts of MHD jets before speciflcally focusing on the 
acceleration (Sec. 3) and collimation (Sec. 4) mechanisms. In Sec. 5, we apply our results to 
GRBs and discuss ultra-relativistic jets within the standard internal-external shock model 
of GRBs. In Sec. 6, we isolate a mechanism for jet instability occurring at large radii (be- 
yond the photosphere) for the case of extreme Lorentz factor and associate it with possible 
dissipation mechanisms. We summarize our results in Sec. 7. 



2. Basic concepts of relativistic MHD jets 
2.1. Magnetohydrodynamic jets 

The scenario of MHD jet formation (Blandford & Payne 1982; Pudritz & Norman 1983; 

Camenzind 1986) can be summarized as follows. The jet is initiated as a slow wind from 
the inner disk by a process which is not yet completely understood, in particular its time- 
dependent character. Most probably, some disk instability is responsible for ejecting the 
matter in the direction perpendicular to the disk surface. The disk wind is flrst launched 
magneto-centrifugally and further accelerated and (self-) coUimated into a narrow beam by 
Lorentz forces. 

Two key parameters which determine the dynamics of relativistic MHD jet flows, (i) 
the plasma magnetization a and (ii) the light cylinder of the magnetosphere i?L- These will 
be discussed below. 

Considerations of stationary MHD flows have revealed that relativistic jets must be 
strongly magnetized (Michel 1969; Camenzind 1986; Li et al. 1992; Fendt & Camenzind 
1996). In that case, the available magnetic energy can be transfered over a small amount 
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of mass with high kinetic energy. On the other hand, a very strong magnetization may be 
in conflict with the MHD assumption lacking a sufficient large amount of electric charges 
which are needed to drive the electric current system (see §3.6). 

Theoretical modelling of jet formation requires to solve the governing MHD equations. 
However, due to the complexity of MHD and the astrophysical boundary conditions indi- 
cated, a completely self-consistent MHD solution for the jet formation process being com- 
patible with all the generic features (MHD self-collimation, accretion-ejection mechanism, 
magnetic field generation, spatial and time scales etc.) does not yet exist. 

In this paper, we will concentrate on solutions to the stationary, axisymmetric, ideal 
MHD equations in the relativistic limit. Then the jet magnetic field distribution can be 
described by a bunch of nested axisymmetric surfaces, measuring the magnetic flux through 
a circular area around the symmetry axis. ^ = (l/27r) J Bp ■ dA. The stationary approach 
has the advantage to obtain a global solution for the relativistic jet MHD structure (e.g. Li 
(1993); Fendt (1997)) on spatial scales and with a resolution which cannot (yet) be reached by 
time-dependent simulations (Kudoh et al. 1998; Koide et al. 2000). Long-term (Newtonian) 
MHD simulations, however, possess the abihty to demonstrate the self-collimation of MHD 
jets (Ouyed & Pudritz 1997; Ouyed et al. 2003). 



2.2. The magnetization pcirameter 

The essential parameter for MHD jets is the magnetization parameter (Michel 1969), 

(7 = ^^. 1 

4Mc3 ^ ' 

The iso-rotation parameter f2F(^) is frequently interpreted as the angular velocity of the 
magnetic field lines. The function $ = Bpr"^ is a measure of the magnetic field distribution 
(see Li (1993)), and M = npVpE? is the mass fiow rate within the fiux surface. Equation (1) 
demonstrates that the launch of a highly relativistic (i.e. highly magnetized) jet essentially 
requires at least one of three conditions - a rapid rotation, a strong magnetic field and/or a 
comparatively low mass load. 

In the case of a spherical outfiow ($ = const) with negligible gas pressure one may 
derive the Michel scaling between the asymptotic Lorentz factor and the flow magnetization. 

Too = <7V3 (2) 

(Michel 1969). Assuming a constant mass flux across a jet with magnetic flux $jet and mass 
flow rate Mjet, the Michel scaling gives V^o — (fip ^)^^^(-^et)~^'^^- It must be mentioned 
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already that Eq. (2) relies on further constraints. For a general relation r^{a) the influence 
of coUimation, gravity and gas pressure^ must be considered. 

Depending on the exact magnetic field distribution $(r, 2;), in a collimating jet the matter 
can be substantially accelerated beyond the fast point magnetosonic point (Begelman & Li 
1994; Fendt & Camenzind 1996), as it is moved from infinity to a finite radius of several 
Alfven radii. As a result, the power law index in Eq. (2) can be different from the Michel- 
scaling (see Sec. 3.5; Fendt & Camenzind (1996); Vlahakis & Konigl (2001)). If the $(r; ^) 
is decreasing outwards, the asymptotic flow is dominated by the kinetic energy (Begelman 
& Li 1994). Still, the speed of the flow at the fast point follows the Michel scaling Eq. (2), 
Ffm = ct(v1>)V3. 

The function <l>(r; z{'$)) = <l>(r; \E') describes the opening of the magnetic flux surfaces 
\E'(r, z) comparable to the action of a "magnetic nozzle" (Camenzind 1989; Li et al. 1992). 
A flux function $ constant along the fleld lines as applied in the Michel scaling corresponds 
to a constant opening angle of the magnetic field. 



2.3. Relativistic features 

At the light cylinder (hereafter I.e.) the velocity of the magnetic fleld lines "rotating" 
with angular velocity f)F(^) coincides with the speed of light. The I.e. is located at the 
cylindrical radius 

7?l(^) = c/npi^). (3) 
For a differential rotation of the fleld line foot points as in an accretion disk, 

the cylinder deforms into a light surface (of a priori unknown shape)^. The I.e. has 
to be interpreted as the Alfven surface in the limit of vanishing matter density (force-free 
limit). Outside the I.e. the magnetic fleld lines "rotate" faster than the speed of light^. The 
existence of a strong toroidal fleld component allows the matter, being frozen into the fleld, 
to slide along the fleld which guarantees v < c also in the region R > Ri^. 

The location of the I.e. determines the relativistic character of the magneto sphere. If 
the light cylinder is comparable to the dimensions of the object investigated, a relativistic 

^In the case of a "hot wind" , the magnetization is not a free parameter anymore. 

^Notc that general relativistic effects affect the shape of the light cylinder. Frame dragging close to a 
rotating black hole implies a second light surface. 

'^As the field line is not a physical object, the laws of physics are not violated. 
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treatment of MHD is required. 

Contrary to Newtonian MHD, in the relativistic case electric fields cannot be neglected. 
The poloidal electric field component is directed perpendicular to the magnetic fiux surface. 
Its strength scales with the I.e. radius, Ep — E± — {r/Ri,)Bp. As a consequence of Ep ~ Bp, 
the effective magnetic pressure can be lowered by a substantial amount (Begelman & Li 
Begelman & Li (1994)). 

A further difference between relativistic and Newtonian MHD is the fact that the 
poloidal Alfven speed ua becomes complex ior r > Rj^, u\ Bp (1 — [r/Ri^Y) — Bp — E\. 
Therefore, Alfven waves cannot propagate beyond the I.e. and only fast magnetosonic waves 
are able to exchange information across the jet. 

Note that the I.e. arises as a relativistic effect due to the rapid rotation of the mag- 
netosphere. This has to be distinguished from the fast proper motion of matter in poloidal 
direction leading to relativistic effects affecting the inertial forces. For relativistic MHD jets, 
both features are interrelated. A rapidly propagating jet must originate in a rapidly rotating 
source. This interrelation is parameterized by the Michel-scaling. 

The I.e. is essentially a special relativistic feature. Close to the black hole horizon 

general relativity becomes relevant. The existence of the I.e. as a natural length scale in 
relativistic MHD is not consistent with the assumption of a self- similar ]ci structure. The 
latter holds even more when general relativistic effects are considered (see Sec. 4). 



2.4. Magnetization of GRB jets 

In this paper we discuss the possibihty that GRBs are generated by ultra-relativistic 
MHD jets. Essentially, two model scenarios for the jet origin may be considered, a jet 
launched from an accretion disk or by a highly magnetized neutron star. As the properties 
of the central source driving the jet arc yet unknown, any estimate on the jet magnetization - 
magnetic field distribution, mass flow rate and rotation - must rely on somewhat hypothetical 
parameters. In the following, we will assume that the jet is launched from the accretion disk 
around a collapsed object and estimate the jet magnetization by constraints from disk theory. 
We do not consider the origin of the disk magnetic field. It could be advected within the 
disk from the ambient medium, or generated by a disk dynamo. 

We constrain the jet magnetization at the foot point radius of the jet a = {Bp^r'^Vlp / c^Mjet) , 
from the disk equipartition field strength B^^ ~ Pgas- This value hmits the toroidal magnetic 
field component which can be built up by the disk differential rotation, B^ < ^eq and also 
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for the poloidal magnetic field amplification by a dynamo process, Bp < Bt- In the case of 
a self-similar advection dominated disk with accretion rate M^cc, we obtain 

-1/2 / M V^/^ 



Seq ~ 7.8 X 10^ G 



Me j \rs 

where Me = 10 Lq/c^ = 2.2 x 1O~^(M/M0) Mgyr"^ is the Eddington accretion rate, a is 
the disk viscosity parameter, and rs the Schwarzschild radius (Narayan et al. (1998), see also 
Fendt & Greiner (2001)). An optically thin standard accretion disk with Thomson opacity 
gives a similar value. A number of GRB jet models discussed in the literature consider 
magnetic field strengths of up to lO^^G (e.g. Lyutikov & Blackman (2001)). Equation (4) 
shows that such a high field strength can never be expected from a disk magnetic field for sub- 
Eddington accretion rates. Hyper-accreting disks around black holes have been discussed by 
Popham et al. (1999). These models provide accretion rates of 0.01 to 10 MqS"^ and consider 
efficient cooling by neutrino losses. As a consequence, the equipartition field strength may 
reach 10^"^ — lO^^G. We note however that a is unknown for hyper-accreting disks. 

The other parameter in the magnetization is the jet mass flow rate. As this is unknown, 
we scale the jet mass flow rate in terms of the disk accretion rate, Mjct — 10~^Macc- An 
upper limit for the jet magnetization can be derived considering the marginally stable orbit 
as jet origin, = rms- Here, we find the largest disk magnetic field strength, together 
with the most rapid rotation, I^f — 1-4 x 10^ (M/M0)"^(r*/3rs)"'^/^ (in the Schwarzschild 
case). For a maximally rotating black hole, (a/M) ~ 1, we have — '^g, where rg is the 
gravitational radius of the black hole. With the maximum jet magnetic fiux constrained 
by the disk equipartition field, B^, ~ B^q, the jet magnetization essentially depends on two 
parameters - the mass ejection rate and the jet origin. 

From the Michel-scaling, Eq. (2), we derive a minimum asymptotic jet Lorentz factor Poo — 
10. 

GRB afterglow observations indicate a total baryonic mass in the burst of about 1O~^M0 
(Piran 1999), implying a hypothetical jet mass flow rate of about Mjet — IO^^MqS^^ if the 
burst lasts for 10 s. Such mass flow rates can never be achieved for disk accretion rates 
constrained by the Eddington limit. 



M,. . 1.4xlO-"MeS-'(^)-'(iL) 



(6) 
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The hyper-accreting stages of accretion disks discussed by Popham et al. (1999) are a way out 
of this dilemma. With Mjet — 10~^Macc, the inferred accretion rate could be < IO^^MqS"^ 
(necessary for neutrino cooling) and the equipartition field strength is increased substantially. 
However, the flow magnetization governing the asymptotic speed will remain about the same. 



3. Acceleration — the asymptotic Lorentz factor 

Assuming axisymmetry and stationarity, the equations of ideal MHD can be re-written 
into two equations describing the force-balance along the field (the MHD wind equation, 
hereafter WE) and the force-balance across the field (the Gmd-Shafmnov equation, hereafter 
GSE). In general, both equations are interrelated as the source term of the GSE depends 
on the dynamics of the MHD wind solution. In turn, the wind acceleration depends on the 
magnetic field structure which is given by the solution of the GSE. However, in the case of 
highly relativistic (i.e. highly magnetized) jets, the infiuence of the moving matter on the 
magnetic field can be neglected and the field structure may be calculated by the force-free 
GSE (see Sec. 4). 

In this section we present ultra-relativistic MHD solutions of the wind equation. For 
simplicity, we will consider the cold wind equation in Minkowski space-time, which gives us 
the freedom to investigate the fiow dynamics for a different choice of magnetization. For a 
hot wind MHD solution in Kerr metric we refer to Fendt & Greiner (2001). 

As the kinetic time scale for the GRB jet propagation, Tkin, is well above the time scale 
when the jet crosses the collimation region, TcoU, 

10i?L 10^\m 
Tcoll ^ ^SX IQ-^S « Tkin = ^ 33 s 7 

c c 
stationarity may indeed be applied for the jet formation region. 



3.1. The relativistic MHD wind equation 

Combining the MHD equation of motion with the conservation laws for energy E, an- 
gular momentum L, magnetization a and iso-rotation Qp, we obtain the wind equation for 
the poloidal velocity tip = Tvp/c in Minkowski space-time, 

4 + 1 = (8) 
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a;2(l - M2 - x2)2 

(Camenzind 1986) . Here, x\ — (f^p L/E) defines the Alfven radius xx and M\ — {A.-K^n'Up) / 
the Alfven Mach number Ma {n' is the proper particle density). In the cold wind hmit the 
wind equation simplifies to a polynomial equation of degree of four for the poloidal velocity. 
The polynomial coefficients explicitly depend on the magnetization a, the flux tube function 
$, and the flow parameters energy E and angular momentum L (see Fendt & Camenzind 
(1996); Fendt & Greiner (2001)). 

At the magnetosonic points the wind equation becomes singular. A finite solution only 
exists if numerator and denominator vanishes together. For this critical wind solution the 
poloidal velocity of the matter equals the speed of the magnetosonic waves at the magne- 
tosonic points. We consider such a solution as a global solution as it is accelerating from low 
velocities at small radii to large speed at large radii. 

The cold wind solution is defined by the following parameter set. The magnetization 
a and the iso-rotation Qp are free parameters and may be constraint by the astrophysical 
boundary conditions. The total energy density of the flow E is constrained by the regularity 
condition at the fast point. With Qp and E also the total angular momentum flow L = E/Qp 
is determined. 



3.2. The ultra-relativistic asymptotic MHD jet 

Figure 1 shows a sample of MHD jet solutions for a parameter set as motivated above^. 
We have calculated the flow dynamics for different magnetization a = 1000 and a = 5000, 
and for a different magnetic field distribution $(r; ^') ~ r^^ with q = 0.01,0.1,0.2. The 
hght cyhnder is located at Rj, — lO^cm. Gravity is unimportant in the cold hmit. Thus, the 
degree of coUimation does not change the character of the solution as a function of radius r 
(e.g. up{r)). However, for a collimated flow, for each radius r the distance from the source 
z{r) is increasing with increasing degree of coUimation. 

Our solutions demonstrate that ultra-relativistic velocities can be achieved by a MHD 



^The figure shows the two positive solution branches out of the set of four numerical solutions of the wind 
equation at each radial point. The branch which starts with low velocity at small radius and accelerates to 
high speed at large radii is the continuous branch of the (stationary) physical wind solution. The other one 
decelerates with radius and is not defined for each radius. For an appropriate parameter set, the latter branch 
may turn into a continuous accretion branch. The intersections of these branches identify the magnetosonic 
points (the Alfven point at i? ~ iiL and fast magnetosonic point at R> Rh) which determine the solution. 
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jet if it is highly magnetized. We obtain velocities up to i^p ~ F ~ 3000 for a — 5000. Prom 
the sequence of plots in Fig.l it can be seen that magnetization a and field distribution 
(parameter q) play an equally important role concerning the jet acceleration. A similar gain 
in asymptotic velocity can be achieved by either increasing the magnetization by a factor of 
five or increasing q from close to the Michel value to 0.1 (compare upper right with middle 
left figure). This confirms earlier results by Begelman & Li (1994) and Fendt & Camenzind 
(1996). 

3.3. Pitch angle of the magnetic field 

It is well known from MHD wind theory that the motion of matter along poloidal 
magnetic field lines also implies induction of a strong toroidal field component - so strong 
that it overcomes the poloidal component for radii larger then the Alfven radius. This has 
been shown early in the case of non-relativistic self-similar MHD disk winds (Blandford & 
Payne 1982), but holds also for jet flows launched from the disk around a rotating black hole 
(Fendt & Greiner 2001). For the MHD jet solutions presented in the present paper we flnd 
the same result (not shown). In particular, the latter paper shows a power law distribution 
for the toroidal magnetic field decay in radial direction for radii larger than the light cylinder, 
i?T ~ The poloidal magnetic field strength, however, decays faster, as g < for any 
reasonable field distribution. For a monopolc type field we have g = and Bp ~ while for 
a dipolar field Bp ~ and q = —1. For the solutions presented here, for the "asymptotic" 
domain (i.e. for r = 10^, z ~ 10^) the numerical solution gives B^/Bp ~ 10^. 

Essentially, this proves again the well known characteristics of a magnetohydrodynamic 
jet which is asymptotically dominated by the toroidal magnetic field component. This may 
have important implications e.g. for the magnetic field distribution in the asymptotic shocks 
and for the interpretation of the polarization structure in the afterglow observations (see 
Sec. 6). 

3.4. Energy balance along the flow 

MHD jets essentially live from the exchange of magnetic and kinetic energy. When the 
jet is launched with low velocity, the energy content is mainly in the magnetic part, i.e. we 
have a Poynting dominated flow. As the flow accelerates, it gains kinetic energy converting 
Poynting flux into kinetic energy flux by Lorentz forces. 

Figure 2 shows the energy partitioning of our MHD jet solutions for high magnetization 
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(7 = 5000. We see that only in the case of ^ ~ the flow remains Poynting dominated also 
for large radii. In the case of a faster magnetic flux divergence {q > 0), the flow accelerates 
substantially beyond the fast surface, converting more and more Poynting flux into kinetic 
energy. For the chosen magnetization the kinetic energy becomes substantial only beyond 
the fast magnetosonic point. Eventually, the energy distribution in the asymptotic flow is in 
rough equipartition between the kinetic and magnetic contributions. This is consistent with 
the claim of Begelman & Li (1994), but also with (Newtonian) numerical simulations of jet 
formation (Ouyed & Pudritz 1997; Fendt & Cemelijic 2002; Ouyed et al. 2003). 

3.5. Asymptotic velocity — a modified Michel scaling 

Another feature demonstrating the influence of a divergent magnetic field is the corre- 
lation (Michel scaling) between asymptotic jet velocity and magnetization (Fig. 3). 

As a measure for the asymptotic {z oo) poloidal velocity "Up^ooi we have taken the 
velocity at the r = 10^ -Rl- The actual asymptotic value might be a factor of two larger, 
however, what is important is that the velocity profile saturates significantly beyond 10 - 100 
I.e. radii. This implies that, as far as the efficiency of magnetic acceleration is concerned, 
we do not have to consider the far distant region of the flow. 

Figure 3 shows three curves. The lowest curve corresponds to the original Michel scaling 
Up,oo ^ cr^^^ in the case of a magnetic field distribution with g ~ 0. If the magnetic fiux 
decreasing faster {q > 0), there is a substantial gain in asymptotic velocity. For the solutions 
with the decreasing magnetic flux (two upper curves in Fig. 3) we flnd a modified Michel 
scaling, comparable to an almost linear relation Up^^ ~ Aa, where A depends on the choice of 
q (this is crucial since large variations of the jet's Lorentz factor are possible in this solution; 
see 6.1). For the model parameters discussed here, we flnd A ~ 10~^/^ for q — 0.1 and 
A ~ 10-^5 for q = 0.2. 

Note that the interrelation Up^^{a) shown in Fig. 3 essentially provides a link between 
an asymptotic, "observable" jet parameter (velocity) and a parameter which is intrinsic to 
the jet origin (magnetization) . 

3.6. Applicability of the MHD approximation 

With the advantage of our approach of knowing the exact solution of the MHD equations 
along the coUimating jet flow, we are able to check self-consistently a posteriori whether the 
approximation of magnetohydrodynamics is satisfled within the calculated flow. The problem 
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is hidden in the fact that for a cold MHD jet one may find arbitrarily high velocities for an 
arbitrarily high flow magnetization. However, an arbitrarily high magnetization may be in 
conflict with the intrinsic MHD condition under which the solution has been calculated and 
which requires a sufficient density of charged particles in order to be able to drive the electric 
current system (Michel 1969). Below a critical particle density the concept of MHD breaks 
down. As a good measure for this critical density we consider the Goldreich- Julian density 
uq (Goldreich & Julian 1969; Lyutikov & Blackman 2001). 

Following the notation of equation (8), the particle density n in terms of the Goldreich- 
Julian density uq as a function of radius along the magnetic field line is 



(see Appendix) . The index ★ denotes a number value at the foot point radius of the magnetic 
field line Vi,. The magnetic field z-component decreases along the opening magnetic flux 
surfaces while the Alfven Mach number increases as the jet accelerates. In the Appendix 
(Fig. 5) we show some example curves of the relative density {n{r)/nQ{r)) derived for the jet 
solutions shown in Fig. 1. Figure 5 demonstrates that the break-down of the MHD concept 
is critically important rather for highly magnetized jets with weak absolute Poynting flux 
(implying a low mass flow rate). For all solutions presented in this paper n/nc stays larger 
than about 1000. 



In this section we discuss the structure of the collimating GRB jet. We apply a stationary 
approach and present axisymmetric solutions of the MHD equations in Kerr metric. Time- 
dependent general relativistic MHD simulations of jet formation in the literature have failed 
so far to span time periods of more than some rotational periods. 

The axisymmetric magnetic field structure of a stationary collimating MHD jet follows 
from the solution of the Grad-Shafranov equation (see above). The problem is to obtain a 
global solution which, at the same time, also considers the local force-balance. So far, fully 
self-consistent solutions of the stationary relativistic MHD equations have not yet been able 
to obtain. The assumption of a self-similar MHD jet is in general a powerful approach in order 




(9) 



4. CoUimation and the compactness problem 



-13- 



to obtain self-consistent MHD jet solutions^ (see Vlahakis & Konigl (2003) for an application 
to GRB). However, for the case of relativistic jets this assumption implies a certain angular 
velocity at the foot point of the field lines fiF(^) ~ (Li ct al. 1992) which is in clear 
contradiction both with the Keplerian rotation of a disk or with the rigid rotation of a central 
body. Similarly, also the magnetization a must be fiux-independent. It is therefore essential 
to treat the relativistic MHD jet in a non self-similar, fully two-dimensional approach. So far 
this has been possible only in the limit of force-free force-balance, neglecting the inertial back- 
reaction of the matter on the field structure (Camenzind 1987; Fendt 1997; Ghosh 2000). It 
is clear that the previous made comments on a self-similar jet structure become even more 
valid in the case of a general relativistic treatment of MHD jets. Resulting solutions for 
"self-similar relativistic MHD jets" ( Contopoulos 1994; Vlahakis & Konigl 2001, 2003) will 
not be free of these problems. 



4.1. Structure of relativistic MHD jets from rotating black holes 

Here we summarize the essential steps in calculating the axisymmetric force-balance of 
relativistic MHD jets (for a detailed discussion see Fendt (1997)). In difference from the 
previous section we now consider the governing MHD equations in the framework of general 
relativity. 

We apply Boyer-Lindquist coordinates in the 3-1-1 split of space-time around a rotating 
black hole of mass M and angular momentum per unit mass, a — J/Mc with the line element 

ds" = a'c^dt' - Qj^ (dcj) - ujdtf - (pVA) dr^ - dO^. (10) 

t denotes a global time in which the system is stationary, is the angle around the axis of 
symmetry, and r, 6 are similar to their fiat space counterpart spherical coordinates^ Here, 
uj is the "frame dragging" angular velocity of an observer with zero angular momentum 
(ZAMO), UJ = {d(j)/dt)zAMO- The lapse function a describes the lapse of the proper time r 
in the ZAMO system to the global time t, a = {dr / dt)zAMO- 



^Still, one has to keep in mind that self-similarity implies further constraints to the solution, (i) It does not 
allow to include the jet axis in the treatment, (ii) Self-similar jets have an infinite radius. As noted already 
by Blandford & Payne (1982) the radially self-similar assumption becomes (iii) "increasingly artificial" when 
the jet has formed at large distances from the disk 

^Thc parameters of the metric tensor are defined as usual, p^= -\- cos^ 0, A= — IGMrjc? -\- a?, 
uj= 2aGMrlcY?, Y?= (r^ + a'^ f - a^Asin^ 6, 
u)= (S/p) sin 6*, a= pTA/S 
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We define the axisymmetric magnetic flux through a loop of the KiUing vector m = 
a;^V0 as 

^ir,e)^^ [ Bp-dA, = X (11) 

(see Sec. 2.1). The indices P and T denote the poloidal and toroidal components of a vector. 
Similar to Eq. (11) the total poloidal electric current is defined by / = — J ajp ■ dA — 

— ♦ ^ -* — ♦ 

— ^aujB^ . In a force-free magnetosphere with = pcE + -j x B the poloidal electric current 
is parallel to the poloidal magnetic field Bp \\ jp and is a conserved quantity along the 
magnetic flux surfaces, I — 

The axisymmetric force-balance perpendicular to the magnetic field is described by the 
Grad-Shafranov equation (hereafter GSE). In the limit of highly relativistic (i.e. highly mag- 
netized) jets, inertial forces have a negligible influence on the structure of the magnetosphere 
and we may apply the force-free hmit of the GSE, 

^V.f«i^i?^V*U-|4^. (12) 

\ j 2aCj d^ ^ ' 

For simplicity, differential rotation of the jet basis has been neglected, I^f(^) = const. = flp. 

~ 1/2 

The two light surfaces are located at the radial position u;l = (±tt/(^F ~ ^^)) • The + sign 
holds for the outer light surface with Qp > uj, while the — sign stands for the inner light 
surface, where Qp < uj. The asymptotic radius of the outer light surface (cZ^l for -2 ~^ 00)1 the 
light cylinder, is denoted by -Rl- Normalizing the GSE (12) using a; — > i?L cD; V ^ (l/i?L) V ; 
^ "if max ^ and / — > /max I is numerically advantageous, but also provides insight in the 
physical characteristics of the solution. The coupling constant gi measures the strength of 
the source term of the GSE, 

c2^Lx " ^lO^^A^ ^WcJ 402iGcm2^ " ^^^^ 

Basically, gi determines the strength of the electric current /max- A high coupling constant 
corresponds to a strong poloidal electric current (respectively a strong toroidal magnetic 
field) and therefore imphes a strong jet coUimation. Note that the coupling between the 
source term and the poloidal field structure considers only electromagnetic quantities. The 
solution is, however, calculated in Kerr metric, and considers also gravity. The link between 
the two governing length scales defined in the problem - the asymptotic light cylinder and 
the gravitational radius - is made by by choosing in terms of the black hole rotation Qh 
(see below). 

Equation (13) apphes parameter estimates as discussed above for GRB jets launched by 
black hole accretion disks. The maximum poloidal electric current in the disk-jet system is 
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given by the disk equipartition field strength. For an advection dominated disk (see Eq. (4)) 
we have S^^'^^ < B^^ ~ lO'^G and I^^^ < 3 ~ lO^^A. Since B^^"" < St ^ Bcq in the 

disk, we derive a maximum magnetic flux from the disk \E'inax ^ (3rg)^i?™^^ ~ lO^^Gcm^. 
The expression for gi can be further simphfled considering that the hght cyhnder radius in 
Eq. (13) is governed by the rotation of the jet foot point Vl-p. As both B^^^ and B^^^ are 
limited by i?eq as a function of radius, we flnally obtain 



For the case of hyper-accreting disks discussed above, the limiting equipartition field strength 
constrains the poloidal and toroidal field equivalently. As a consequence, the coupling gi 
remains the same. 



Here we present an example solution for the stationary axisymmetric force-free mag- 
netic field structure of GRB jets calculated from the Grad-Shafranov equation (12). The 
GSE is solved applying the method of finite elements. This allows to consider the complex 
geometrical structure of the black hole - disk - jet system. Details of the method of solution 
are discussed elsewhere (Fendt 1997). We summarize the following characteristics of the 
solution. 

- The shape of the collimating outer jet boundary is not known a priori. It is determined 
by the internal structure of the jet and is particularly constrained by the regularity 
condition along the outer light surface (i.e. a smooth transition of the magnetic field 
lines) . 

- The "free function" /(^) in the GSE source term is chosen from the analytical solution 
of an asymptotic jet in perfect coUimation. 

- The interrelation between the two characteristic length scales of the solution, the 
asymptotic light cylinder and the gravitational radius - hence, the link between elec- 
trodynamics and gravity - is made by choosing the angular rotation Of(^)- Here, 
Vl-p = OAQn corresponding to the Keplerian rotation at the marginally stable orbit. 
With that Rl = lOrg for a = 0.8. 

- Other parameters are the magnetic fiux from the black hole in terms of the total fiux, 
^Bn/^max = 0.2, and the profile of the disk magnetic fiux distribution. 




(14) 



4.2. Example solution for a fully collimated GRB jet 
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Figure 4 shows the resulting structure of the coUimating magnetic flux surfaces. The solution 
extends from the inner light surface to the asymptotic regime of a fully coUimated jet. We 
find a perfect collimation of the asymptotic jet in agreement with analytical models of (Appl 
& Camenzind 1993). We find a rapid collimation of the outfiow into a cylindrical shape 
within 50rg distance from the central black hole. The expansion rate (asymptotic jet radius 
in respect to the jet foot point radius) is about 10. The asymptotic jet radius is 5 I.e. radii. 
Due to the numerical resolution we do not obtain solutions with larger jet radius. 

Whether the jet is self-collimated by magnetic tension or pressure coUimated by an 
ambient medium cannot be answered by this approach. 

An interesting fact is that the flux surfaces emerging from the black hole (i.e. from 
the inner light surface in our approach) asymptotically partly intersect with the outer light 
surface. This implies that the mass load in the jet provided by the accretion disk is predom- 
inantly in the outer layers of the jet at radii r > i?L- Here, we expect the high velocities 
calculated in the previous section together with a strong toroidal magnetic held. The inner 
part of the jet with flux surfaces emerging from the black hole, should be dilute of matter, 
i.e. we have a hollow jet stream. This part of the jet may be powered by the electromagnetic 
interaction with the black hole itself (Blandford-Znajek mechanism) . 

The strength of the coupling is similar to the case of active galactic nuclei (AGN). As 
the coupling constant Eq. (13) measures the strength of the coUimating toroidal magnetic 
held strength in terms of the de-coUimating effects poloidal held and rotation^. Therefore, 
we may expect a similar structure of the jet magnetosphere in both examples. A strong 
coupling implies efficient collimation; we find that collimation angles as low as 0" — 5" can 
be achieved with proper tuning of parameters. From AGN jets we know that they are well 
coUimated with asymptotic jet radii up to of about 100 i?L- Just from the similarity of the 
coupling constant we might expect the same scaling also for jets in GRB's. To calculate such 
large jet radii is, however, beyond the capability of our numerical code. 

5. The "bottle-neck" instability for highly magnetized jets 

In the previous sections we have presented stationary solutions for highly magnetized 
MHD jets - the fiow dynamics along the coUimating magnetic field and the axisymmetric 
magnetic field structure. Our time scale estimates (Eq. 7) indicate that the stationary ap- 
proach may indeed be applied when considering the jet formation region. Nevertheless, it is 



'Note that the couphng constant corresponds to the inverse magnetization for force-free jets 
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clear that the GRB itself must be the consequence of a time-dependent process. The current 
picture of GRB invokes highly relativistic shells moving with different speed, catching up 
with each other by forming a highly energetic shock front - the gamma ray burst (Piran 
1999). 

Prom our MHD solutions we have found indication of a possible excitation mechanism for 
a flow instability. This instability essentially relies on (i) a high magnetization and (ii) a re- 
collimation of the governing magnetic flux tube. The solutions discussed in Sec. 3 have been 
calculated for a decreasing function $ along the jet corresponding to a increasing opening 
of the flux tube. However, the two-dimensional jet magnetic fleld structure as presented 
in Sec. 4 shows an interesting property. The global jet coUimates into a cylindrical shape 
with the single flux surfaces turning from an initially conical outflow close to the disk into 
an alignment parallel to the jet axis. This implies a re-collimation of the single flux tube - 
easily visible when comparing the width of the flow channel along the outer flux surfaces (see 
also Pendt & Camenzind (1996)). The maximum width is typically located in the region of 
the strongest curvature at a distance of several i?L from the central source. 

Solving the wind equation by considering a re-collimation of the flux tube we flnd 
that no stationary wind solutions exists beyond the radius where re-collimation happens. 
We cannot directly answer the question what exactly happens there as this is beyond our 
stationary approach. Por comparison, one may think of the earth magnetosphere reflecting 
solar wind particles from the polar region or the more general picture of a "magnetic bottle" 
of converging magnetic field lines. 

In the appendix we show a typical MHD wind solution as discussed in Sec. 3, but mod- 
ified for the re-collimation feature by prescribing an increasing magnetic fiux tube function 
$ beyond a certain radius. In this example we applied a cosine-like recollimation of the fiux 
tube $(r) ~ -^i^j-^jj^ certain radius range beyond r > 50rA — 50-Rl (see Eq. Bl). 

Pigure 6 shows an example velocity distribution along the field line for different magneti- 
zation. Por high magnetization the magnetic fleld is strong enough to brake the matter 
sufficiently. The exact dynamical behavior depends on the parameters of the wind solution. 
Essentially, the re-collimation leads to a deceleration of the wind flow before the stationary 
solution ceases to exist. Applying a spatially limited region of re-collimation, e.g. a sinu- 
soidal variation of \l'(r) between two radii, there may exist a unique stationary flow solution 
which can be decelerated below the magnetosonic velocity. However, in reality such a unique 
combination of fleld distribution and mass flow rate may not exist. 

Therefore, we propose the above mentioned indication of a "bottle-neck" instability as 
a generic reason for the velocity variation and the formation of shock waves in GRB. This 
instability is working predominantly in the outer layers of the jet and in the coUimation 
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regime at distances of several I.e. radii from the central object. 

The onset of the instability depends on the details of the flow propagation. However, 
in general, the essential conditions are strong magnetization and re-collimation. Due to the 
Michel scaling, this corresponds in particular to high Lorentz factors as present in jets of 
GRBs. In our example model we compare magnetizations of cr = 10, 1000, 5000 and a certain 

sinusoidal variation of the flux tube function (see Appendix). It can be seen that for low 
magnetization a stationary solution exists from the jet foot point to the asymptotic regime. 
For higher magnetization, the stationary solution ceases to exists at a certain radius. For the 
parameter choice applied here, wc do not find a stationary flow domain around the location 
of local maximum re-collimation. In the case of further recoUimation, no stationary wind 
solution will be found beyond this radius. In our case, with an again increasing size of the 
flux tube, stationary solution can be found beyond the non-stationary regime. This, however, 
has no consequence as the jet flow has become non-stationary already. We note that it is 
not straightforward to estimate the exact location of the radius where the instability occurs. 
However, its seems to occur at large enough radii to assume that it is beyond the Compton 
photosphere. 



6. Ultra-relativistic MHD jets in the GRB picture 
6.1. Varying Lorentz factor and internal shocks 

Mass-load 

Variations in the Lorentz factor of 10 — 1000 (10^ < F < 10^) as suggested by the 
GRB internal shock model are hard to obtain from the original Michel scaling. In fact, 

■ —1/3 

'Up,oo ~ ^jct would imply a variation in the mass load by several orders of magnitude 
more. The existence of a modified Michel scaling as described in §3.5 (see Fig. 3) leaves the 
possibility of having a substantial change on the Lorentz factor by the variation of the initial 
flow magnetization. A variation in the jet magnetization can be caused by a change in the 
mass injection rate from the accretion disk into a jet with temporarily constant magnetic 
field implying Up^^ ~ M^^ . It is interesting to note that these variations are directly induced 
by conditions at the source pending contamination by the surrounding material as discussed 
next. 

Ambient mass-entrainment 

One would expect parts of the jet to be contaminated by the ambient material (Daigne 
& Mochkovitch 2000) as to cause local variation to the intrinsic Lorentz factors (as given 



-19- 



by the disk mass loading or magnetization). In fact, it is likely that the variations in the 
Lorentz factors are induced by a combination of varying mass-load and mass-entrainmcnt. 
This would result in a multiple shock mergers inducing GRBs (as in the internal shock 
scenario) which can be seen if mass entrainment occurs mainly beyond the compactness 
radius. 



6.2. Magnetic energy dissipation 

It is argued that, for models of internal shocks in GRBs to successfully reproduce the 
GRB temporal features, different shells of matter should have a comparable energy and 
their different Lorentz factors should arise due to modulation of the accelerated mass (Piran 
(2001) and references therein). In MHD jets, the energy (mainly magnetic at the base of 
the flow) is roughly constant as the magnetization is not expected to vary much. A variable 
mass-load combined with the entrained mass, as we have said, offers the modulation needed 
to account for the wide range in the Lorentz factors^. By itself, such a variable source is not 
enough to explain the variable light curves since magnetic energy dissipation mechanisms 
will not necessarily be efficient or adequate in reproducing the GRB temporal structure^*^. 
For the Crab Pulsar wind in which energy is transported predominantly as Poynting flux, 
the fluctuating component of the magnetic field in such a flow can in principle be dissipated 
by magnetic reconnection and used to accelerate the flow (Kirk & Skja^raascn 2003). In 
our model, such mechanisms might be at play as a result of the "bottle-neck" instability. A 
conversion of the fraction of the energy of the accelerated particles into radiation is interesting 
within the GRB context. The pulsation of the emitted radiation in our picture, as we have 
argued above, would be linked to the variation in the accretion rate at the source. We note 
however that it is not yet clear how much of the magnetically dissipated energy can be seen 
in the form of 7-rays (Spruit et al. 2001). 



6.3. Collimation and polarization 

The results from §4 suggest that ultra-relativistic MHD jets are highly coUimated. Ex- 
trapolating from what we know in the case of AGN jets this implies that collimation angles 



^This may also lead to a narrow range in the Lorentz factor distribution as part of the jet with higher 
Lorentz values would in principle induce higher entrainment; a notion which remains to be confirmed. 

^"^Dissipation in shocks (internal and external shocks) is fundamentally different from magnetic energy 
dissipation in MHD outflows where non-thermal electron distributions is not guaranteed. 
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as small as 0° — 5" are feasible. Such an efficiency in coUimation is a statement that the jet is 
asymptotically dominated by the toroidal magnetic field component (§3.3). This may have 
important implications e.g. for the magnetic field distribution in the asymptotic shocks and 
for the interpretation of the polarization structure in the afterglow observations. 

7. Summary 

We studied ultra-relativistic MHD jets in the context of GRBs. We have presented 
stationary solutions for the governing MHD equations considering the axisymmetric structure 
of the coUimating jet and the acceleration of matter by the magnetic field. Special and general 
relativistic effects have been taken into account. 

Lorentz factors up to 3000 can be obtained for suffient strong jet magnetization. The 
advantage of our approach of actually solving the governing MHD equations is that we can 
prove a posteriori the applicability of MHD itself by comparing the matter density with the 
Goldrcich- Julian density along the jet flow. The asymptotic jet magnetic fleld is dominated 
by the toroidal component. The energy distribution within the asymptotic jet is almost in 
equipartition between magnetic and kinetic energy. 

For the structure of the axisymmetric MHD jet we flnd a rapid coUimation to almost 
perfect coUimation within a distance from the black hole of about half the jet radius (of 5 
light cylinder radii in this example). 

Among the features that are crucial to standard models of GRBs, are (i) the modifled 
Michel scaling {up oc ~ ^jct^) which allows for a plausible variation in the Lorentz factor 
by variations in the source parameters such as the disk magnetization and/or mass-loading; 
(ii) the high degree of coUimation [6^0° — 5") which is within the range of values derived 
from breaks of afterglow light curves. We isolated a jet instability that develops at extreme 
Lorentz factors. This instability which seems to arise at large radii, if it occurs beyond the 
compactness radius, would result on an optically thin emission once its associated magnetic 
field is dissipated. We expect different magnetic field configurations to lead to different 
dynamics (e.g. asymptotic profiles, dissipation efficiency). This, in principle, could reproduce 
the large diversity and duration range of GRBs. 

Among the open questions are the extreme magnetic fields required, and the effect of 
baryon contamination (via mass-entrainment) on the distribution of the Lorentz factors. It 
is also not clear at this stage why and how the engine (the underlying hole-disk system here) 
should stop for several dozens seconds before bursting again as seen in GRBs. 
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A. Applicability of the MHD concept 

As discussed above highly relativistic jets must be strongly magnetized (Michel (1969); 
Camenzind (1986); Fendt & Camenzind (1996); see Sec. 2. 4). Magnetohydrodynamic jets 
live from magnetic to kinetic energy conversion. More kinetic energy per fluid element can 
be delivered for jets with correspondingly strong magnetic field, or low mass flow density, 
respectively. On the other hand, the MHD concept itself requires a minimum density of 
charged particles in order to allow for electric currents in the flow. The jet flow mass 
density naturally decreases along the flow as the jet originates in a very small region close 
to the central source and first expands in almost radial direction before it collimates into 
a narrow beam. When the flow density is below the critical density at the MHD radius 
'"MHD, the MHD assumption breaks down. With our model approach, we may determine 
self-consistently the possible breakdown of the MHD conception as our MHD wind solution 
delivers all the dynamical and electromagnetic properties of the flow. 

Melatos & Melrose (1996) quantify the range of validity for ideal MHD, E + v x B = 
0, by inspecting the one-fluid equations for a cold, neutral plasma. They consider the 
generalized Ohm's law for fluctuations of time scales r on characteristic length scales A in 
the magnetohydrodynamic quantities. In the end, there are four constraints for ideal MHD. 
First, it is required that j « nev and PcV « j with the charge density Pc- If fhe conduction 
current if dominating the displacement current. Maxwell's equations imply that for MHD 
(Anjr/E) ~ {cr/Xf. With the relation 

-max 1,— ,— «1 (Al) 



ujqX^ \ ' TV ' A ^ 

(with ujq = A/47rne^/mp) Ohm's law actually reduces to ^ + x i? ~ 0. In the case of a 
pulsar wind with A ^ -Rlj uj c:^ Q^, = Qp and r ~ u"^ = Rh/c, Eq. (Al) takes the form 
Fu;^ << Uq which is equivalent to the estimate of Michel (1969). For the Crab pulsar the 
MHD radius is located at tmho — 10^i?L (Melatos & Melrose 1996). 

A similar criterion for the applicability of MHD is a matter density above the Goldreich- 
Julian charge density (Goldreich & Julian 1969; Michel 1969). The Goldreich- Julian particle 
density along a magnetic flux surface ^'(r, z) is 

_ (V^-V)(r^QF) _ g,(r;^) 

and is not negligible for rapid (relativistically) rotating magnetospheres. Again, a particle 
density below the Goldreich- Julian density indicates that not enough charges are available to 
carry the electric current. From Eq. (A2) we see that if the particle density profile decreases 
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f aster than the component, there may exist a critical distance tmrd where the matter 
particle becomes lower than the Goldreich- Julian density. 

Usov (1994) has apphed this approach in order to model GRBs generated by millisecond 
magnetars of 3 x lO^^G dipolar field strength, 

In the case of a leptonic plasma wind Usov derives a critical distance Tmhd — lO^^cm for 
the given stellar magnetic field strength and field distribution and the rotation rate {Qp = 
lO^s""^). It is, however, another question to extrapolate these variables over more than six 
orders of magnitude in distance from the star to the region where the critical radius is 
located in this model. A similar approach has been applied in the context of "Poynting- 
dominated" outflows in GRBs (Lyutikov & Blackman 2001), assuming a certain magnetic 
fleld and density distribution. In this model, large-amplitude electromagnetic waves break 
at the MHD radius located at about lO^^cm while accelerating particles to ultra-relativistic 
speed. 

In difference to the previous work, in our paper we calculate the solution of the rela- 
tivistic MHD wind equation from the jet basis to the asymptotic regime providing us with a 
set of dynamical properties of the flow. In particular, this allows us to check self-consistently 
(but a posteriori) the consistency of our solutions with the MHD conception. With that we 
may also compare different jet geometries - rapid or weak collimation - and may constrain 
the parameters of the jet mass loading and the magnetic field distribution. We derive an 
expression for the matter particle density from the definition of the Alfven Mach number. 



p 



Ml = (A4) 



with the proper particle density n', the gas entropy fj, and the poloidal velocity Up = Tvp/c 
(Camenzind 1987; Fendt & Camenzind 1996). Re-writing Eq. (A4) gives the proper particle 
density along the magnetic field line as a function of the Alfven Mach number, 

n' = -^^^— (A5) 

Relating that to the Goldreich- Julian density Eq. (A2), we find Eq. (9). As the Alfven Mach 
number generally increases faster than the z-component of the magnetic field decreases along 
the flow, the ratio (n'/ncj) will decrease. In Fig. 5 we show the proflle of the density ratio for 
four of the example wind solutions presented in Fig. 1. For all four cases the mass flow density 
stays above the Goldreich- Julian density. The MHD jet parameters are the magnetization 
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(7 = 1000, 5000 and the foot point magnetic field strengtfi of Bp — lO^G, lO^^G. In particular, 
we find that the density ratio approaches a constant value n' /uqj > 100 for large radii. This 
reflects the fact of a jet collirnafAon. For a further (substantial) expansion of the jet flow for 
these choice of initial parameters, the mass density would become under critical. 

B. Recollimation — the "bottle neck" instability 

Here we show example solutions of the MHD wind equation for the case of a local re- 
collimation flux tube function $(r). For this purpose we modified the decreasing magnetic 
flux tube function <l>(r) by simply adding a cosine shape increase within a certain radius 
range along the jet. For the solutions presented in Fig. 6 we have prescribed 

$(r) ~ r~°-^ for r < 50 ta 

$(r) ~ r-°-^(^l + 1.25(^l-cos(^7r))) (Bl) 

for 50rA < r < 2050rA 
$(r) ~ r-°-2 (1 + 1.25 (1 - cos tt)) for r > 2050 ta 

We have computed the wind solution for several different magnetizations and otherwise 
equal parameters. Figure 6 shows that for low magnetization {a — 10) the stationary wind 
solution reaches from the foot point into the asymptotic domain. For higher magnetization 
((T = 1000, 5000) no stationary solution can be found along a certain radial range around the 
point of (local) maximum re-coUimation (see $(r)-plot). 
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Fig. 1. — Relativistic MHD wind/jet solution. Poloidal velocity Up along the magnetic field 
line. Parameter: a = 1000 (left), a = 5000 (right), Parameter: $(r) ~ r'l, q = 0.01 (top), 
g = 0.1 {middle), q — 0.2 (bottom). 
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Fig. 2. — Relativistic MHD wind/jet solution. Magnetic energy in terms of total energy 
as a function of radius along the magnetic field line. Note the fast magnetosonic point as 
intersection of both curves. Parameter: a — 5000, $(r) ~ r~^, q — 0.01, 0.1, 0.2 (from top to 
bottom) . 
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Fig. 3. — Relativistic MHD jet solution. Modified Michel scaling Up^^{a) for a different 
choice of the magnetic field distribution ^(r;^) ~ with q — 0.2 (top curve), q — 0.1 
{middle curve), q — 0.01 {bottom curve). 
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Fig. 4. — The axisymmetric magnetic field structure of a jet from a rotating black hole as 
solution of the Grad-Shafranov equation in Kerr metric ^(r, 2;). The jet originates within 
5rg of the accretion disk. Parameters: BH angular momentum parameter a — 0.8, coupling 
constant gr — 0.717, jet radius rjet = 50 r^, asymptotic light cylinder i?L = 10 r^, black hole 
magnetic flux ^bh = O.l^totai- By prescribing gi, rjet and the disk flux distribution, the 
shape of the coUimating jet boundary {thick curvejis determined by the regularity condition 
at the light surface. The inner boundary is the inner light surface around the ergosphere 
(white hemisphere). 
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Fig. 5. — Relativistic MHD wind/jet solution. Particle density in terms of Goldreich Julian 
particle density as a function of radius along the magnetic field line. Parameter: a = 1000 
(left), a = 5000 (nght), initial magnetic field strength 5^,* = lO^G (top), B^^^ = lO^^G 
( bottom) . 
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Fig. 6. — Relativistic MHD wind/jet solution. Profiles of the magnetic flux tube function 
{top) and the poloidal velocity {below) along the field line for a re-coUimating magnetic flux 
tube. Magnetization a — 10, 5000 (from top to bottom). 



